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A MEMOIR ON THE PROBLEM OF DISTURBED ELLIPTIC 
MOTION. 


[From the Memoirs of the Royal Astronomical Society, vol. XXVIL, 1859, pp. 1—29. 
Read March 9, 1858.] 


I VENTURE to take up the problem of disturbed elliptic motion, for the sake of 
a further elaboration of the analytical theory. The points which present difficulty are 
the measurement of longitudes in the varying plane of the orbit, and (in the lunar 
theory) the determination of the position of the orbit by reference to the varying 
plane of the sun’s orbit; it is, in memoirs and works on the lunar and planetary 
theories, often difficult to discover where or how (or whether at all) account is taken of 
these variations, and the analytical mode of treatment is for the most part very imperfect. 
I must except always Hansens Fundamenta Nova [investigations orbite vere quam 
Luna perlustrat, Gotha 1838] where the points referred to are treated in a perfectly 
rigorous manner. There is, however, a want of clearness in the form under which his 
investigations are presented; and the comprehension of them is greatly facilitated by 
Jacobi’s remarks, published under the title “Auszug zweier Schreiben des Prof. Jacobi 
an Herrn Director Hansen" (Crelle, t. XLII. pp. .12—81 (1851). Jacobi observes 
that the integration of Hansens system of differential equations introduces seven 
arbitrary constants, which, in the expressions for the coordinates referred to fixed axes, 
reduce themselves to six. The seventh constant, neglecting the disturbing forces, is 
in fact a constant which determines the position in the orbit of the arbitrary origin 
from which the longitudes in orbit are reckoned. I have, in my paper “On Hansen's 
Lunar Theory,” Quarterly Mathematical Journal, vol. 1. pp. 112—125 (1855), [163], 
termed this origin “the departure-point, and longitudes measured from it “departures.” 
The seventh constant may be taken to be the departure of the node. I reproduce in 
the present memoir the explanation of what is meant by the departure when the plane 
of the orbit is variable. If the problem is treated by the method of the variation of the 
elements, the seventh constant becomes, like the other elements, variable; and we have 
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thus a seventh variable element, the departure of the node. 'The element just referred 
to (the departure of the node) forms, with the longitude of the node and the inclination, 
a group of three elements, which determine the position of the orbit and of the 
departure-point. 'The coordinates of the planet are in the first instance taken to be 
the radius vector, longitude, and latitude; but the before-mentioned three elements 
being considered as given, the position of the planet depends only on the radius vector 
and the departure. These may be then expressed in terms of the remaining four 
elements; as to the choice of these four elements, it is to be remarked that there is 
one element which only enters through the mean anomaly, and that there is great con- 
venience in representing with Hansen the mean anomaly by a single letter; and that in 
the various formule we may use, in the place of the element implicitly involved in the 
mean anomaly, the mean anomaly itself, or treat the mean anomaly as an element; the 
four elements may be taken to be the semi-axis major, the eccentricity, the mean 
anomaly, and the departure of the pericentre. And joining to these the before-mentioned 
three elements, we have the system of elements represented in the memoir by 
a, € 9, w, c, 0, d It has been assumed so far that the three elements determine 
the position of the orbit and departure-point in reference to a fixed plane and origin of 
longitudes; but we may suppose more generally that, instead of the fixed plane and 
origin of longitudes, we have a variable plane or orbit of reference and a departure- 
point in this variable orbit of reference. The quantities which determine the orbit of 
reference and departure-point are naturally taken to be the departure of the node, 
longitude of the node, and inclination; these are assumed to be given functions of the 
time, and they are in the memoir represented by oa’, 6’, $'. The three elements of the 
planet's orbit (viz. departure of node, longitude of node, and inclination) in relation to 
the orbit of reference and departure-point therein, are in the memoir represented by 
£, ©, ®, and the system of elements ultimately adopted is therefore a, e, g, w, X, O, ®. 
I obtain formule for the variations of these elements under two different modes of 
expression of the disturbing function: first, when the disturbing function is expressed 
in terms of the radius vector and departure and of the three elements X, ©, ®: 
secondly, when the disturbing function is expressed in terms of the seven elements 
a, e, 9, w, 2, O, c. The establishment of the two sets of formule just referred to 
constitutes the chief object of the memoir; but the memoir contains some other inves- 
tigations and formule in relation to the general subject. 


The coordinates of the planet are 
r, the radius vector, 
v, the longitude, 
" the latitude. 


The attractive force at distance unity is for convenience represented by n?a*, which 
denotes, therefore, an absolute constant; but the significations of n and a are not yet 
defined. 


The disturbing function, as used by Lagrange, is denoted by Q, that is Q — — R, 
if R be the disturbing function of the Mécanique Céleste. 
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The equations of motion are 


d dr o ($)- ay na’ dO 
didt | 99 9 Adi r( = 


dt okie i 
7° cos? A do) 
a ? di ~ dy’ 
d 2 : dv? _ dQ 
al ae +r cos ysin y (F) = 


where Q is regarded as a function of r, v, y, or (as this may be expressed) where 
OSO (r, v, y). 


If we neglect the disturbing forces, the planet moves in an ellipse; and taking 
a to represent the semi-axis major, the mean motion will be n. The mean anomaly, 
which I call g, will be a function of the form nt+c; but as c only enters through g, 
it will be convenient to use the mean anomaly g (considered as implicitly involving 
an arbitrary constant c) in the place of an element, and I write 


a, the semi-axis major, 

e, the eccentricity, 

g, the mean anomaly, 

0, the longitude of node, 

¢, the inclination, 

©, the distance of pericentre from node. 


I assume also 
, the true anomaly, 


z, the distance of planet from node, 


S 


, the reduced distance from node. 


We have then r and f given functions of ¢ and the elements, viz. we may write 


r=a elqr (e, 9), \ 


f= elta (e, g), 
(read  elqr. elliptic quotient radius, and elta. elliptic anomaly). These values 
j _ a(l—é) ; 
satisfy r— reme fF Moreover z, z, y, are the hypothenuse, base, and perpendicular 


of a right-angled spherical triangle, the base angle whereof is $; the equations which 
connect these quantities are therefore 

tan z = tan z cos ¢, 

sin y = sin z sin 4, 

tan y = sin 2 tan $, 


COS 2 = COS T COS Y, 
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equivalent, of course, to two equations. The first and second of them give in fact 
æ, y, in terms of z and ¢. 
The value of z is 
z=0 +f, 


so that « and y are given functions, and the longitude v is given in terms of æ and 


0 by the equation 
v=2+8, 


and consequently the three coordinates r, v, y, are by the system of equations given 
in terms of ¢ and the elements. 


From the equations which connect z, æ, y, $, treating all these quantities as 


variable we deduce 
sec? z dx = cos $ sec? z dz — tan z sin $ d$, 


cos y dy —sin $ cos zdz + sin z cos $ d$, 
sec? y dy = tan $ cos «dx + sin «sec $ dd, 
sin zdz =cos y sin æ dæ + cos zsin y dy, 


equivalent of course to two equations; and the system is easily reduced to the more 
convenient form ; 


dz = cos $ sec? y dz — tan z cos’ æ sin $ dẹ, 

dy — sin $ cos c dz + cos g tan z cos ddd, 

dx = cot$ sec g se? y dy — tan z cosec $ d, 

dz — cos ġ dae 4- cos æ sin $ dd, 
joining to these equations the 

dz — dO -- df, 

dv =de + dé, 


and considering at present the mere analytical forms, first if d$ =0, dO —0, we have 
dæ = cos $ sec? y dz, 
dy = sin $ cos «dz, 
dz — df, 
dv =de. 

Next, if dy — 0, dv — 0, we have 

dæ = — tan zcosec $ d, 
dz =-tanz cot $d, 
dz = cos $ da, 
da = — dð, 
dz + cos dé = 0. 


I remark also that the equation, tan y — sin z tan $, may be written in the form 


cos? h sec? y + sin? $ cos? g= 1. 


C FDE: 35 
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The equations 
r=a elqr (e, g), 
f^ elta(e, 9) 


treating all the quantities as variable, give 


die ENT IR Lo da — a cos f de, 


J1—-e 1+ecosf 
quae m dy d SETA tecosf gy 
(1 — ey l-e 
to which is to be joined 
TE ac (l-e) sinf g l-e a(2e-1+ecosf) yy 


EET V * 14.505 ^^ vi (1c ecosf 


all which formule will be useful. 


[212 


If we treat the elements as constant, then in the foregoing expressions for dr 
and df, we must attend only to the part involving dg, and must put this equal to 
ndt; the values first obtained for dæ, dy, dz, dv, correspond to this assumption, and 


we have 
pa _ nae sin f 
ET RC 
wer 
dt Pray’ 
dz | na V1—é 
a r? ; 
da cos ġ sec? n ER 
dt ? - $ y T? ? - 
dv _ 3 na? V1 — 
d; 7909 $ sec p : 
dy _ naæNl -e 
gg in $ cos " , 
and we then deduce 
d dr _ nate cos if 
dt dt es Pn 
3i (r costy St) = 0, 
2,44 An M 
: G A = — cos? d sin y sec? y M 1 
NOI S EER I S 
dt \dt E 73 à 


values which satisfy the undisturbed equations. 
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The disturbed equations may be dealt with in the usual manner by the method 


of the variation of the elements, and attending only to the variations of the elements 
we have 


dr =0, 
dv — 0, 
dy — 0, 
dr dQ 
dw dQ 
d (r® cost y Ej we A dt, 
dy\ dQ, 
a(n ài) = ay & 
or, what is the same thing, 
dr =0, 
dv =0, 
dy =0, 
naesinf dQ 
PERO dr s 


dna 16 cos $ — 92. dt, 


do 


dna? 1— esin $ cos s = du 


dt, 


where as before Q=Q(r, v, y). 
In virtue of the relations dv — 0, dy — 0, we have the above-mentioned equations, 


da = — tan z cosec $ d, 
dz — —tanzcot dd, 


dz=  cosó da, 
da — — d6, 
dr +cos¢dé=0, 
we have 
d sin $ cos æ = — sin sina da + cos æ cos $ dẹ, 
= cos gcos $ sec? z d, 
= secacos psec’ y dó ; 


and the last two equations for the variations become 
dna? /1— e cos $ — na? /1— e sin ddd = —- di, 


Je D dQ 
dno? V1 — @ sin $ cos z -- na? N1 — e! sec æ cos $ sec? y d = et 


35—2 
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and attending to the equations cos? d sec? y + sin? ġ cos 5-1 we deduce at once 
g q y 


T— ,, 20 at dQ 
dna?Vl—e- cos $ sec y. deqonndousm gm 
do -— (— sin $ cost a SP dt + cos ġ cosa dt). 


Now the position of the planet may be determined by the quantities v, z, 0, h, or 
we may consider Q as a function of the last-mentioned quantities And if on the 
right-hand side Q =Q (r, v, y) as before, the formule of transformation are 


I. +o av +T dy= 2 ar +O de +O +55 do 
where 
dv = cos ¢ sec? y dz — tan z cos? æ sin $ dp + dé, 
dy =sin $ cos «dz + tan z cos «cos $ dẹ, 
and we have 


dQ dO 

dr dr’ 

dO do 

dð dv’ 

dQ, ,.dQ dQ 
ag 7 ian (—sin g eos py; C008 cosa iy) 
dQ „dQ dO 
a ( cos $ sec yrs 7 + 8in densa), 


where on the left-hand side Q =Q (r, z, 0, $); and these equations give 


dQ dQ 
cot 2 Ta = cot b gy — cosee $ ga $ 


an equation which is satisfied by O =Q (r, z, 0, $). We have thus \ 


dr=0 
dv = 0, 
dy = 0, 
naesinf _ dQ 
i i dr % 
dnav1—@ = a 
dr 
cot z dQ 


SS | wi di, 
na? V1—e do 
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which may be replaced by 


dr'= 0, 

do = © cosec $ - do di 
nN- e dé ^" 
_ —cobt$ dQ di 
naàjNl—e d$ ^ 

gei _ dQ 
VI -e dr 7^ 
dna VIZE = 20 gy 

dz 


dp m cot z dQ 
di na? V1— e d$ " 


where as before Q=0 (r, z, 0, $). 


I remark that in the case of any central force whatever, we have an element h 
corresponding to na? V1— e in the elliptic theory, and the system for the variations is 


dr — 9, 

qp = 2560 # dé di, 
dz = ie ag d^ 
gu oe 
d$ = ut T dt, 


where Q=Q (r, 2, 0, $). 
Imagine a point in the orbit, which I call the departure-point, the angular 
distances from this point are termed departures. And I write 
p, the departure of planet, 
w, the departure of pericentre, 
c, the departure of node, 
pot 


£ —-bp-oc, 


so that we have 


C=a-<a. 


www.rcin.org.pl 


278 A MEMOIR ON THE PROBLEM OF DISTURBED ELLIPTIC MOTION. [212 


I write also 


5, the longitude in orbit of departure-point, or, as it may be termed, the adjustment; 


that is 


8-0—oc. 


In the undisturbed motion the departure-point is simply a fixed point in the 
orbit, but when the orbit is variable, the departure-point is taken to be the point 
of intersection of the orbit with any orthogonal trajectory of the successive positions 
of the orbit, a definition which is expressed analytically by the equation, 


do = cos $ dé. 
The equation, z 2p — c, gives 
dz = dp — de = dp — cos $ d6, 


or, what is the same thing, 


dp — dz -- cos $ dé. 


But we have dz+cos¢d@=0, and consequently dp —0, an equation which expresses 
that the increment of departure, in so far as such increment arises from the variation 
of the elements, is equal to zero. Or, what is the same thing, the total increment 
of departure is equal to the infinitesimal angle between two consecutive radius vectors 
of the planet. \ 


I propose to consider the departure-point as a point which is constantly defined 
as above, viz, when the orbit is variable, the departure-point is the point of inter- 
section of the orbit with any orthogonal trajectory of the successive positions of the 
orbit; and as a particular case of the definition, when the orbit is- fixed, the 
departure-point is simply a fixed point on the orbit. The orbit here considered is 
that of the planet and the position of the planet is determined by the departure and 
radius vector (the latitude being zero), and this is assumed to be the case whenever 
the departure is spoken of, and it is such departure which is denoted by the letter |. 
But we might consider a departure-point (defined as above) upon any other orbit 
whatever, and use such departure-point as an origin of longitude (for instance, in the 
lunar theory we might consider a longitude measured along the variable plane of the 
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sun’s orbit from a departure-point, defined as above, in that orbit) and the position 
of the planet would then be determined by means of the longitude, latitude, and 
radius vector. The term sidereal longitude is, I think, used in Physical Astronomy 
rather loosely to denote the longitude in the mean ecliptic from the mean equinox, 
less the precession; so defined it is not practically different from, and may I think 
in all cases be replaced by the longitude as measured from a departure-point in the 
mean ecliptic. 


Returning from this digression, the assumed equation, de= cos d, gives the 
expression for the variation de of the departure of the node, and we now have in 
the place of the former six equations the seven equations 


dr — 0, 
dp — 0, 
resin _ 10 v 
v1l—@ dr 
GE dQ 
dnaæNli-e= d; d^ 
tz dQ 
FIND LU NA A 
? naæNl -e dọ 
td do 
de SE dt, 
na l—e do 
dO, 
dp e QU FALL 
na? 1—e dó 


where as before Q =Q (r, z, 6, $). 


But the value of z is z=p—o, and Q can be expressed, and that in a single 
way only, viz. by means of the substitution of p—o in the place of z, in the form 
Q, — Q (r, p, c, 0, $), and if on the right-hand side €) — OQ (r, z, 0, $) as before, then we 
have 


io. do 
dr dr’ 
io do 
E EM ae 
.dQ  dQ 
de (dz 
io do 
dé dd’ 
dO dO 
dO  d0" 
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where on the left-hand side Q =Q (r, p, c, 0, $) The function Q so expressed 
satisfies, of course, the partial differential equation 


(which conversely implies that p, e only enters through the function p— 6), and it 
also satisfies the partial differential equation obtained from the before-mentioned equation 


dQ dQ dQ 
ots =cot $ -Jz — cosec ó dà ' (Q — Q (r, z, 0, $)), 
by the introduction of the transformed expressions of the differential coefficients, and 


which may be written 
dQ 


£y dQ dQ 
o =— cot T — cosec $ d8^ 


where Q =O (r, b, c, 0, $). 


Using the last-mentioned equation to transform the value of dd, the expressions 
for the variations become 


dr = 0, 
dp — 0, 
nae sin f _ dQ dt 
/1—e de TA 
dna?N1—6e un dt, 
y 
- “do cosecd dQ 
dip = 4 TLI: EE TNI EN ied 
? ng? 1— e do na? 1—6g dé a 
d cot $ do i 
2 najNl—e dp ’ 
cosecd dQ 


dé dt, 


s na? V1—e dd 
where Q=Q (r, p, c, 0, p) as before. 


I suppose now that the orbit of the planet, instead of being referred to a fixed 
plane, is referred to a moveable plane or orbit of reference. It is assumed that the 
longitudes in the orbit of reference are measured from a departure-point defined as 
above,—that is, from the point in which the orbit of reference is intersected by any 
orthogonal trajectory of the successive positions of the orbit of reference. And the 
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position in regard to the fixed plane, of the orbit of reference, and of the departure- 


point in this orbit, are determined by 0, o’, $'—that is, we have for the orbit of 
reference, 


0', the longitude of node, 
c', the departure of node, 


¢’, the inclination. 


The position of the planet’s orbit in relation to the moveable orbit of refer- 
ence is determined in like manner by 9, >, ®,—that is, we have for the planet’s orbit 
in relation to the orbit of reference, 


©, the longitude of node, 


X, the departure of node, 


®, the inclination. 


Hence if, as before, 0, c, $, belong to the orbit of the planet considered in relation 
to the fixed plane, 2—o, 0 —6', 0—0, will be the sides of a spherical triangle, the 
opposite angles of which are $', 180° — $ and d. 


Putting for shortness S- È ee S' 20 —c', G —0 —0', so that these symbols denote 
S, the distance of node, along planet's orbit, from fixed plane, 
S', the distance of node, along orbit of reference, from fixed plane, 
G, the distance in fixed plane of the nodes on fixed plane, 
the sides of the spherical triangle are S, S’, G, and the opposite angles are ¢’, 180^ — $, ®. 
OQ. HI. 36 
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Calling the sides A, B, C, and the opposite angles a, b, c, the general formule for a 
spherical triangle give without difficulty, 


dC=  -—cosbdA — cos a. qp ., 2n 4 sin B sino ; 


sin C 
sin b eos C sin a sin A cos b 
dud O alee Cord DO ip ctii 
eos C sin b sin @ sin B cos a 
db = sin C dA sin C dB sin C ; 
and conversely 
dc = cos B da + cos A 45,4 80.5 T bein ag 
sin C 
FT sin B ja cos c sin A 3 dá sina cos B ig 
sin c sin c sin c 
dB = — ®t B Ja sin Á gb _ nin bos A 25 
sin c sin c sin c 
which, in the present case, become 
+ 2 , zı, , Singsin ¢’ sin G 
do = cos S do cos S' dg’ + "une" CI ae dG , 
uis sin S cos $ sin S' ,,, sin $ cos S 
BeN eurer cue fc 
_ €08 $ sin S end eee sin $' cos S' 
piis sin $ mt ar T sin $ ud 
and 
ee , 7 sin S sin Ss’ sin d 
dG = cos $ ds’ — cos $ dS + peu uL. 
en sin $' cos G sin $ sin S cos ¢’ 
xil v dua cuu RESI C M 
,_ _ cosGsing’ zy sin $ sin S cos $ ] 
dp =- mne at "NAM T Er CS Wo 
and we have also 
dS =d% —doe, 
dS’ —d6 — dc', 
dG —d0 —d6'. 


Hence, observing that sinSsin$ = sinS'sin$/, the preceding equations may be 
written 
dd = cos S d$ — cos S'd$'-- sin S sin $ d0 — sin S'sin Q' dO 
dO = dc'— cosec Psin Sd -- cot $ sin S' dd’ + cosec ® sin $' cos S (d0 — dé’), 
d£ — de — cot sin S d$ + cosec $ sin S' d$ + cosec ® sin ¢’ cos Y (db — de’), 
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and 
d@ =d0 + cos $' (dO — do’) — cos $ (dX — do) + sin S sing d®, 
dp = -—cosec G sin ¢' (d@—do’)+ cot G sin $ (dX — do) + cosec G sin S' cos $' db, 
dó = -— cot G sin $' (dO — do^) + cosec G sin ġ (dX — de) + cosec G sin S cos dẹ, 


and it is proper to remark, that in obtaining these equations no use has been made 
of the equations de — cos $ d0, do’ — cos $' dé’. 


The term in dO which contains do’, &c. may be written 
(do’ — cos $' d") + cos $' d& + cot «b sin S' do’ — cosec ® sin $ cos S dé’, 
which is equal to 
(do — cos $' d0") + cot ® (sin S dd’ — cos S' sin $' dé’); 
and the term in d£ which contains de, &c. may be written 
(do — cos $ d0) + cos $ d0 — cot P sin S d + cosec ® sin $' dd, 
which is equal to 
(do — cos $ d0) — cot ® (sin S d$ — cos S sin $ d0); 
reductions which depend on 
cos $' — cosec ® sin $ cos S = — cot ® sin $' cos S', 
cos @ +cosec ® sin $' cos S' = — cot sin $ cos S, 
or, what is the same thing, 
cos $' sin ® — sin $'cos cos S = sin $ cos S, 
cos $. sin  — sin $ cos ® cos S = — sin $' cos S", 


which are relations between the sides and angles of the spherical triangle. And we 
then have 


d= (cos S d$ + sin S sin $ dé) — (cos S' d -- sin S' sin $' d0’), 
d =(do’ — cos ¢’d6’) — cosec ® (sin S d — cos S sin $ d0) -- cot ® (sin S' d$ — cos S' sin $' dé’), 
dX-(de —cos$d0)-— cot 4b (sin Sdo — cos S sin $ d0) + cosec ® (sin S' d$' — cos S' sin $' d0’), 
expressions which may be simplified by omitting the terms (de — cos ¢’ dé’) and 
(de — cos ġ d0) — 

Next substituting for do, d$, d0, their values, we obtain, 


dr = 0, 

dp =0, 
maesinf _ a0 y 
"E Mr 


36—2 
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denied alii d dt, 


dO dO. | Dagmar- . ww dh 
Md cun — ie (sins Bag ne (cot $ -Jg tcosec ó A) dt— (sin S’d¢’+sinS’sin ¢’ dé’), 
cot ® dO gan) a dek nå E r 
gr one (s3 co je S (cot ó dg ^ Cose $ 2) dt-+cosec (sin S'd$' —cosS'sin $' dé’), 
do 
de Seid 3 s (eese B sin S (cot $ exi cosec 4) dt+ cotd(sinS'd$' —cosS'sin $'d6), 


where Q — Q (r, p, c, 0, $), as before. 

But Q may be expressed in the form Q=Q (r, y, X, ©, 6, c', 0', $), or dis- 
regarding o’, 0', ¢’, in the form Q=Q (r, p, =, ©, 4), and to effect the transformation 
of the differential coefficients we must write, 

d® — cos S d$ + sin S sin $ dé, 

de = ` — eosec ® (sin S d — cos S sin ¢ d0), 

dX, = (de — cos $ d0) — cot ® (sin S d$ — cos S sin $ dé), 
or, what is the same thing, 

d$ = cos S do — sin S sin ® dO, 

dé = cosec ó (sin S d® + cos S sin & dO), 

do = dX, — cos d + cot $ (sin S d® + cos S sin  d9), 


and substituting in 


ae dO, 
; dr tratta à I + Ta dé 
S dO, 
ap + = 2 d ub =o P 5 de + 16 d, 
if on the right-hand side Q= O (r, p, e, 0, p) as before, then we have 
io do 
UP" wn 
dQ dQ 
dp dp’ 
dQ dQ 
dX da’ 
^e C cos ® + cot $ cos S sin o) d s + cosec $ cos Ssin 77 -sinSsnd T. 
TS = in SS 4 cosec ¢ sin S cn at cos S T , 


where on the left-hand side Q = Q (r, p, E, ©, D). 
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The last three equations give 


Eu o mm dQ dON dQ 
Seger + sin S (cot $ + cosee $ 35) = jb 
. gad dQ dQ, dQ dQ 
d" tk cos 8 (cot $ aor cosec $ T) =cot DoS + cosec 565, 
and the formule for the variations become 
dr S 
dp = 0, 
naesinf dQ, 
g Vl- ar x: 


e Ls 
dna? V1—e EE dt, 


—cot® dO cosee 5 dQ 


do S OTTERSU NES dt — (cos S’d¢’ + sin S' sin $'d6^), 
dX = LL oe dt +cosec ® (sin S'do’ — cos S' sin gd), 
do = E = dt + cot ® (sin S'do’ — cos S' sin $'d6^), 


where Q 2 Q (r, p, 2, ©, o). It will be recollected that the value of S’ is =@-— ø. 
It may be noticed that 


dX, — cos ® dO = sin ® (sin S' dd’ — cos S' sin g’ dé’). 


The system just obtained is, except as regards the terms involving d$' and dé, 
precisely similar in its form to that in which the planet is referred to a fixed plane, 
or where Q =O (r, p, c, 0, $), and this is of course as it should be. 


We have now 
p- 5A 
r-aelqr (e, g), 
f= elta (e, 9), 
so that the position of the planet is determined by means of the elements a, e, g, w, X, 9, ®. 
To find the variations of these elements, substituting for r its value in terms of f, 
the first, third, and fourth equations are 


all-e) _ 
1+ecosf ’ 
naesinf — dQ 
Vv1-é dr 


dnev1—e= ES dt, 
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which give 


(day Madea aa Bs = n conf g QE CE PHI d Fa 
= 2 p, 
e cos f df -- sin fde m Lus EN Retr gine a 
l-e /1—e dO 
di de 4- ys da — nate dp dt, 
and we thence obtain 
2e sin f dQ 2(1+ecosf) dQ 
da = ——— —— dt ———— a dt, 
3 nV1—é dr “+ na(1—&) dp 
de — vi € inf IRAE de, 
na na? 1—e6e dp 
|^ N1— e cos f dO _ (+ecosf) sin f dO 
fe nae aR Tne a ERA TU OR 


the last of which equations, combined with 


dra (ke oon y ia Qua 


Q-ey de 
gives 
Paler du hacer) ee d _  Qrecosf)sinf dQ g 
nae (1 4- e cos f^) nae . ri 


The fourth equation of the formule for the variations, viz, dp — 0, gives 0 — da + df, 
and therefore de — — df, that is, 


/1—ecos f dQ di 


nae dr 


(2+ecosf)sin f dQ 


"ur naàWl—e dp 


za di, 


and the complete system becomes therefore 


2esinf dQ 2(1-F-ecosfy dQ 
da = dt + dt, 
nV1—@ dr.  aa(1—e) dp 
oe V1—-@sinf dO y, e+2cosf+ecos'f dO, 
na dr na? V1— e dp 
— (1-4) (<2 + cos f+ ecos'f) dO y_ (2-ecosf)sinf dO y 
nae (1 4- e cos f) dr nare dp 
P is poe (2 +ecos/) sinf dQ s, 
nae dr nae dp 
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—coto dO cosec® dQ 


EN. as. "ad ua (cos S" do’ + sin S' sin o d^, 
cot b dQ : " , ps ! JA’ 

ae qe + cosec ® (sin S’ dd’ — cos S' sin $' dé’), 
cosecD dO - Qual ARUM DT 

da- UNI OA dE == di + cot (sin S' dø’ — coss’ sin ¢’ dé’), 


where, as before, Q = Q (r, p, Z, ©, $). This is the first form of the expressions for 
the variations of the elements. 


But we may in the disturbing function Q replace r, p by their values in terms 
of a, e, g, w, and if on the right-hand side Q has the last preceding value, we have 


dQ dQ dQ, +22 dQ dQ dQ 


qs dat Gg de Ut Te MEN 7o 10 -- 12 do 
-2 dQ dQ 
d; 9" mE LIE dz 3 49 + 73 d, 


where on the left-hand side Q =Q (a, e, g, w, X, ©, ©); and the expressions . for. the 
differentials dr and dp, are 


UAE Aoi a ae sin f 
gi It ta a cos f de + js dg, 
n$ [ke con pain f ree 

and we have therefore 
do — 1-e do 
da = 1+ecosf dr’ 
aD, — dQ „@+ecosf)sinf dQ 
dé Tui A ge l-e dp ’ 
dQ _ c (1 -- e cos f£ dQ 
ug c MIS d* u—ef dp’ 
dO _ do 
de dp’ 
aD, _ dO, 
do | dD’ 
dQ, _ dQ 
dE dz,’ 
dO _ dO, 
d® dO’ 


where on the left-hand side Q = Q (a, e, g, w, È, ©, 9), 
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. and we thence obtain for the variations the new system of formule, 


2 dQ à 
"URL 
1-e dO Vv1-e dO 
nae dg needa ’ 
2 dO. 1-e dQ 
T ey na?e de 7^ 
Vv1—e dO 
nate de ' 


— cot dO cosec ® dO 


s dt 


nevi-ed= ^ mgwWi-ede 
2$ (cos S'do’ -- sin S' sin $'d0^, 
cot dQ 
IRSE; ur dd 
+ cosec (sin S'do’ — cos S' sin $'d6^), 


cosec® dQ 
na»N1 —e d^ 


+ cot ® (sin S'do’ — cos S' sin $'d6^), 


where, as before, Q=Q (a, e, g, w, Z, ©, d). This is the second form of the expressions 
for the variations of the elements. It is hardly necessary to remark, that if in either 
system of formule we omit the terms involving d@’ and dọ’, and in the place of 
Z, ©, ®, write c, 0, gd, we have the formule for the variation of the elements when 
the orbit of the planet is referred to a fixed plane, and the disturbing function is 
given under the form Q=Q (r, p, o, 0, $), or Q=Q (a, e, 9, w, c, 0, $). 


The demonstration of the two preceding systems forms, as before remarked, the 
object of the present Memoir. But it is proper to give also the systems for the 
variations of the elements in the form in which they would have been obtained, if 
the notion of the departure had not been introduced into the investigation. To do this 
I revert to à preceding system of equations, which may be written 


dr z, 
d went = dt, 
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br teh dO 
dnai- e= d; 9^ 
-— dQ coto dQ cotz dQ 
d Pe cosec @ d y exi a 
á naV1—e dÓ na Nl- e dz (= na? VI — e dd ) 


where Q2 Q (r, 2, 0, $). 


Substituting in these equations for v, z, the values aue) and O--f, we 


1+ecosf 
obtain 
2esinf dQ 2 (1-- ecos £y dQ 
da = = ——di+ ——_—__> 4- dt, 
= mired f ma-e dr 
~l — esin f dQ e + 2 cos f + e cos? f dQ 
fuw c m s aie, dé ^ 
d (1 — &) (— 2e + cos f+ e cos? f) dQ 3, (2 +e cos f) sin f dQ y 
g nae (1 4- e cos f) dr na?e te 
V1 — e! cos f dO (2 + e cos f) sin f dO | coto dO 
© = — ——— —— di + —.di— —— 
à nae dr 4 d na? V1— e dz na? l—e dd a 
cotó dQ cosecd dO 
= — MR —— ad Te oe er 
a na2Vi—e dz na? V1 — e: d0 p 
dQ 
Bp cosec ġ © d 
nat Vi-@ d$ ^^ 


where Q=Q (1, z, 0, d), as before; and which is the first system for the variations. 
of the six elements, a, e, g, ©, 0, >. 

But if in the disturbing function we replace r, z by their values, then if on 
the left-hand side Q=Q (r, z, 0, $), as before, we find 


1-é do eR 
1+ecosf dr zog 
= 4, Qr e cos f) sinf dQ dQ 

ERE UE l-e dz de’ 
ae sin f " (1-ecosfy dO dO 
Vica dr TM Qs X dz dg’ 
dQ _ dO, 
do | d$' 

dQ, _ dQ 
dz ~ dG’ 

dQ _ dO, 
dé ~ ab’ 


where on the right-hand side Q=Q (a, e, g, ©, 0, $). 
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The formule for the variations thus become 


da Am 
na dg 
m 1-edO , Wi-edo 
«flats nae dg nde dO ” 
2 dO 1—edOQ 
No uet RS 
Vi-edQ , cote dO 
E "nde de “ “navi adg © 
dO 
TT WEE DAE ELEM 
? na? V1 — e dO 
im us cosec $ dQ y 
naNvi—e@dp ^" 


_ _cosecd dQ 
na? V1 — e d0 


[212 


dt, 


. where, as before, Q=Q (a, e, g, ©, 0, $). This is the second system of formule for 


the variations of the six elements a, e, g, ©, 0, ¢. 


The last-mentioned system may be easily deduced from Jacobi's canonical system 


of formule, viz. putting 


9[, the constant of vis viva, 

B, the constant of areas, 

G, the constant of the reduced area, 

F, the constant attached to the time, 

(9, the angular distance of pericentre from node, 


H, the longitude of node; 


then the canonical system is 


dX = ay d 
dB- fpi 
dG = a5 d^ 
d =- Sy di, 
46 - — 72 at 
d$ - — 20 at 
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and the expressions for 2, &c. in terms of the elements a, e, g, ©, 0, $, are 
Y = — iw, 
B= na?V1—e, 
C= na? V1—e cos 4, 


1 
= a Sarg” 


bx HM 
- u 


and the transformation can be effected without the slightest difficulty. 


I shall conclude with the demonstration of a formula which occurs implicitly in 
Hansen's Lunar Theory, and which may probably be useful for other purposes. 


Write 
p, the radius vector, 7 for t, 
a, the true anomaly, r for t, 
that is, let p, 4, be what r, f, become when the time f£, in so far as it enters 
explicitly in g, and not through the variable elements, is replaced by an arbitrary 
quantity, 7. 
And suppose, in like manner, 
y, the mean anomaly, r for t, 
A, the departure, r for t, 


and let lp denote the logarithm of p; then we have, attending only to the variation 
of t, in so far as it enters through the variable elements, 


e: (2 + e cos f) sin f (1 4- e cos f£? 
gom NIS TT, CT nt IE Ld 
Pe (2 + e cos ^) sin y (1 +e cos y} 
dy = tg eae a Uem ^ 


al (1 + e cos q^) cos yr (1 +e cos y) e sin yr 
dlp = a da — Pane HXERAT tos de re ee ae E dy. 


We hence deduce, 


dip: E520 T. (ar dar) 


a (1— e) 
uds, Sen ein f (3 -becon ) — oper — Be Poo ra e sin y» (1 4- e cos f£ 
a (1 — e) (1 +e cos ^) (l—e)F(1+ecosp) "' 
37—2 
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the coefficient of dy being zero. And substituting for da, de, dg their values, viz. 


2esin f dQ 2(1+ecos fy dO 
= - —— dt + -- di, 
ne nV1—e dr H na(l—ey dp 
/1-— e sin f dO e + 2 cos f + e cos? f dQ 
2 Nead nn [We d a 000 gor J 0 
» «^ du de abe Te 
(1 — &) (— 2e + cos f+ e cos? f) 20 qp 4 (2 +ecos f) sinf dO gy 
nae (1 + ecos/) nave d 


the equation becomes 


dlp + PS E (df — dy) 


-V1-e . dQ 
EET T ee " 
1 


+ pape +e cos yr — cos (f — y) 2+ econ /)| Te dt 


But we have p=a+f, A=a+%, and, consequently, f— y — p — 4, and the equation 
becomes 


dip + an $ (dp — dà) 


—VI-¢ 
nai seeps 95 T 


1 


ie pps aa N) cos) di, 


dp 


which is the equation referred to; the expression on the right-hand side, omitting 
the factor dt, is, in fact, the portion not involving the arbitrary functions II, I', of 
Hansen’s function R (Fund. p. 43), viz. it is in Hansen’s notation, 


an dQ 
|B cos (o,— 2) - 1+ zee aa ww [eos - —X)— 1 qu We 


am , i9 
T " dr? 


= sin (y, X) 7—— 


where a, n, e, r, p, v, ^à, Q (Hansen), correspond to a, n, e, r, p, p, ^, n? a^* Q, of the 
present Memoir. 


2, Stone Buildings, W. C., 3 March, 1858. 
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